Spectral properties of a resonator driven by a superconducting single-electron 

transistor 
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We analyze the spectral properties of a resonator coupled to a superconducting single electron 
transistor (SSET) close to the Josephson quasiparticle resonance. Focussing on the regime where 
the resonator is driven into a limit-cycle state by the SSET, we investigate the behavior of the 
resonator linewidth and the energy relaxation rate which control the widths of the main features in 
the resonator spectra. We find that the linewidth becomes very narrow in the limit-cycle regime, 
where it is dominated by a slow phase diffusion process, as in a laser. The overall phase diffusion rate 
is determined by a combination of direct phase diffusion and the effect of amplitude fluctuations 
which affect the phase because the resonator frequency is amplitude dependent. For sufficiently 
strong couplings we find that a regime emerges where the phase diffusion is no longer minimized 
when the average resonator energy is maximized. Finally we show that the current noise of the 
SSET provides a way of measuring both the linewidth and energy relaxation rate. 



I. INTRODUCTION 

When a resonator is coupled to a superconducting 
single-electron transistor (SSET) tuned close to the 
Josephson quasiparticle (JQP) resonance the flow of 
charges can either damp the resonator motion, poten- 
tially leading to coolingfi*^ or pump i^dn^i^d^ leading 
to laser-like states of self-sustaining oscillation, depend- 
ing on the choice of operating point. Recent experiments 
using a superconducting stripline resonator coupled ca- 
pacitively to the SSET island demonstrated a laser ef- 
fectj^ whilst the opposite effect of cooling was demon- 
strated using a nanomechanical resonator Similar cool- 
ing effects and laser like instabilities occur close to other 
transport resonances in the SSET-resonator systemfii^iiii 
as well as in apparently rather different systems such as 
a driven optical cavity coupled parametrically to a me- 
chanical resonatorfi^ii^iii Furthermore, useful analogies 
can be made^ '^^i^^ with quantum optical systems such as 
the micromaser where a cavity resonator interacts with 
a sequence of two level atoms 

The SSET is tuned to the JQP resonance by appro- 
priate choices of the drain source and gate voltages ap- 
plie d^^i^^'^" (see Fig. [T] for a schematic illustration of the 
JQP cycle). Close to the resonance the charge dynamics 
of the SSET island is similar to a driven two-level system 
coupled to a bathj ^^i^^ Charge is transported through 
the system via a combination of coherent tunneling of a 
Cooper-pair and two successive quasiparticle decay pro- 
cesses. The center of the resonance occurs when two 
states of the SSET island differing by a single Cooper- 
pair, |0) and |2) have the same charging energy. For op- 
erating points where the electrostatic energy of state |2) 
is less than that of the |0) state the SSET tends to emit 
energy to the resonatorJi^ The decay processes in the 
SSET generate a current whose average value and fluc- 
tuations provide a natural source of information about 
the dynamics of the resonator i^^ii 

For sufficiently strong coupling, the energy emitted to 
the resonator can lead to a variety of different limit- 
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FIG. f; (a) Superconducting single electron transistor: two 
tunnel junctions and a gate capacitor Cg form the SSET is- 
land to which the resonator is coupled capacitively. A drain 
source voltage, Vds, and gate voltage, V^, are used to tune the 
operating point of the device, (b) JQP cycle; Cooper-pair 
tunneling at the left hand junction between island states |0) 
and 1 2) is interrupted by two quasiparticle tunneling events 
which take the island charge back to |0) via the |f) state. 



cycle statesjii^ii^ The existence of the limit-cycle states 
is shown clearly in the steady state properties of the res- 
onator's density matrix. However, information about the 
important dynamical time-scales of the system such as 
the resonator linewidth and energy relaxation rate is only 
obtained by going beyond the steady-state of the system 
to examine the spectrum of fluctuations that occur about 
this state. 

In this paper we use a combination of numerical and 
analytical methods to investigate the spectral properties 
of a resonator pumped by a SSET tuned close to the 
JQP resonance. Direct numerical evaluation of the rel- 
evant spectra allows us to obtain both the energy re- 
laxation rate and linewidth of the resonator. We find 
that except within transition regions, the peaks in the 
resonator spectra have a Lorentzian shape with widths 
that correspond to the real parts of particular individual 
eigenvalues of the system. For weak to moderate cou- 
plings the resonator linewidth behaves in a way that is 
typical of self-sustained oscillators)^ below threshold it 
is simply half the energy relaxation rate of the resonator, 
while above threshold, in the limit-cycle regime, it be- 
comes much narrower. We show that the linewidth in 
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the limit-cycle regime is set by the phase diffusion rate 
of the system in a way that is similar to a laseri^ For 
sufficiently strong coupling we find that a more complex 
behavior emerges and the operating points (of the SSET) 
at which the linewidth is narrowest no longer correspond 
to the point where the resonator energy is maximized, 
an effect which was also seen in a recent study of a very 
similar systemji We further show that the full current 
noise spectrum of the SSET gives direct access to the en- 
ergy relaxation rate^i and, within a limit-cycle state, the 
phase diffusion rate of the resonator. 

This work is organized as follows. Sec.|TT]is devoted to 
the simple model we use to describe the SSET-resonator 
system and the methods used in its solution. We begin by 
outlining the Born-Markov master equation for the sys- 
tem, in which the evolution of the reduced density matrix 
is given by a Liouvillian superoperator, and then sum- 
marize the steady state properties of the resonator. Next 
we describe the method used to calculate the resonator 
spectra numerically. In Sec. IIIII we give the numerically 
calculated fluctuation spectra of the resonator. Following 
a characterization of the features observed we go on to 
show how the eigenvalues of the Liouvillian matrix can 
be used to obtain the widths of the peaks observed in 
the spectra. The calculation of the eigenvalues is signifi- 
cantly less challenging numerically than a full calculation 
of the spectra, allowing us to then investigate the behav- 
ior of the linewidth for a range of coupling strengths. 
To understand the observed behavior of the linewidth we 
make use of an analytical approximation in Sec. IIVI and 
thereby link the eigenvalues to the relevant physical pro- 
cesses in the system. Finally, in Sec. |V] we show that 
the current noise spectrum of the system can be used to 
measure the linewidth and energy relaxation rate of the 
resonator. The Appendix contains additional details of 
the analytical calculations. 



II. MODEL 

A. Master Equation 

The Born-Markov master equation describing the evo- 
lution of the reduced density matrix, p{t), of the SSET 
and resonator in the vicinity of the JQP resonance is 
given hyr^ 

p{t) = [i^ccPW] + -CqpPlt) + C^pit) 

= Cp{t). (1) 

The first term describes the coherent evolution of the 
density matrix under the Hamiltonian iJco, 

i/eo - Ai?|2)(2|-:|^(|0)(2| + |2)(0|) + |^ 

+ ^mn'^x^ + mn'^Xsx{\l){l\ + 2|2)(2|) , (2) 



where AE' is the detuning from the JQP resonance and 
Ej is the Josephson coupling energy. The resonator 
has frequency fi, effective mass m, momentum operator 
p, position operator x and Xg parameterizes the SSET- 
resonator coupling strength. The second and third terms 
of Eq. [T] describe the dissipative effects of quasiparticle 
tunneling and the resonator's environment respectively, 

£qpP(t) = r|o)(i|p(i)|i)(o| + r|i)(2|p(t)|2)(i| 

-^{|l)(l| + |2)(2|,p(0}, (3) 
■Cdp(i) = ^^^^^ ^ 1^ [x,p{t)]] 

(4) 

where F is the quasiparticle tunneling rate and {• , •} is 
the anti-commutator. The effects of the resonator's en- 
vironment are parameterized by the external damping 
rate 7cxt and the thermal occupation number next = 
^gSO/feBTcxt _ with Text the temperature. We have 

neglected the (weak) dependence of F on the position of 
the resonator and the difference in the tunneling rates of 
the two quasiparticle decay processesJ^ 

For the numerical analysis of the system we use a Liou- 
ville space representatiou j^^i^^i^^i^^'^^ following the nota- 
tion introduced in Ref. 0. The Liouvillian C, appear- 
ing in Eq. [1] can be expressed in terms of an eigenvector 
expansion, 

oo 

where \rp)) are the right hand eigenvectors, {{lp\ the left 
hand eigenvectors and Xp the associated eigenvalues. The 
steady-state density matrix in Hilbert space pss is equiva- 
lent to the right hand eigenvector in Liouville space cor- 
responding to the eigenvalue Aq = 0, i.e. pss ^ ko))- 
Multiplication on the left by {{Iq \ is the equivalent of the 
Hilbert space trace operation)^ {{^o\p{t))) = Tr[p(t)]. 

An appropriate truncation of the oscillator basis allows 
Eq. ^ to be solved numerically to find the relevant eigen- 
values and eigenvectors of £. We do this by using the 
Matlab implementation of the ARPACK linear solver 
In order to use a large number of resonator states we 
neglect the parts of the density matrix corresponding to 
coherences involving the |1) state since these are decou- 
pled from the charge states of the system Additionally 
we make the approximation that the coherence between 
resonator energy levels with a large separation in energy 
can be neglectedi^i^ 

B. Steady state behavior 

We now briefly review the steady-state properties of 
the resonator as a function of the dimensionless coupling 

strength k = ^° and the detuning AE as this will 
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provide important points of reference for the study of 
the spectral properties in the following sections. For con- 
creteness we focus our analysis on the regime of a high 
frequency resonator in comparison to the relaxation time 
of the SSET (Jl ^ F), which is achievable in experiments 
using a superconducting stripline resonator.— This com- 
plements previous work on the low frequency regimcji*^ 
although we note that our methods are also suitable for 
a slow resonator. For this kind of resonator the tem- 
perature can be sufficiently low for thermal effects to be 
unimportant and so we take next = throughout. We 
also choose to work in a regime where the Josephson 
coupling is relatively weak compared to the quasiparticle 
decay rate and hence we take F — Vds/eRj. We use a 
value for the Josephson energy of Ej = 1/16 eVds, with 
a junction resistance Rj ~ h/e". Working in this regime 
has the advantage that it is also possible to develop an 
approximate analytical description of the dynamics f^^i^^ 
Furthermore, when the quasiparticle decay rate is rela- 
tively large this should be the dominant source of dephas- 
ing for the SSET island charge and so we do not need to 
consider additional environmental effects. 

When ^ F the interaction between the SSET and 
resonator is quite weak except for near values of /S.E 
where there is a matching of the eigenenergy of the SSET 
to the energy level separation in the resonator, 
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khVL = ±Jae^ 



(6) 



where fc is a non-zero integer and the sign should match 
the sign of AE. For AE < (fc < 0) these resonances 
correspond to the transfer of energy from the SSET to 
the resonator. For sufhcient coupling the resonator is 
driven into states of self-sustained oscillationsi^ 

The steady state of the resonator is nicely described 
by the distribution P{n) = Tr[|n)(n|/5ss], where \n) is a 
Fock state of the resonator. The resonator is said to be 
in a fixed point state when the P{n) distribution has a 
single peak at n = 0. We define a well-defined limit- 
cycle state as corresponding to a P^n) distribution with 
a peak at 7^ and additionally a small value for P{0) 
(for the purposes of plotting we choose 'small' to mean 
-P(O) < 1 X 10~'^). Finally we define the transition re- 
gion between the two states. In the regimes studied here 
this occurs via a continuous transition, starting from a 
fixed point state increasing the coupling k leads first to 
a wider P{n) distribution and then the peak in the dis- 
tribution moves away from n = 0. We define the tran- 
sition region as when the peak in P{n) is at n 7^ but 
there is also substantial weight at P(0) (again wc choose 
P{0) > 1 X 10~^). Although this particular definition 
of the transition region is of course somewhat arbitrary, 
it nevertheless proves a useful indicator in what follows. 
The resonator can display bistable (and multistable) be- 
havior in this system)^ but here we focus only on a pa- 
rameter range where this is not the case and the Pin) 
distribution only ever has a single peak. 

The steady state properties of the resonator, (n) = 
(a^a), where a is the resonator lowering operator, and 
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FIG. 2: (Color online) (a) Average resonator energy (n) and 
(b) Fano factor Fn, for varying AE and n, with = 10 and 
7cxt = 3 X 10~* (we adopt units such that F = 1). Within 
the inner dashed line the resonator is in a limit-cycle state, 
between the two lines it is in the transition region (as defined 
in the text) and elsewhere it is in a fixed point state. 



Fn ~ ip^) I ('''■) I with n = n — (71), are shown in Figs. [2^ 
and[2b for varying detunings and coupling strengths. The 
plots arc centered around the k = —1 resonance [Eq. ©], 
which for our parameters occurs at AE — — l.SOeVds- 
Similar plots were given in Ref. 0, but here we go to 
higher coupling. As the coupling is increased we see that 
(n) increases on resonance up to a maximum value for 
K ~ 0.005 and then decreases. Fn shows a large peak at 
the point where the peak in the P(n) distribution moves 
from n = to n, 7^ [see Fig. The value of F„ 

drops as the coupling is increased further and the system 
develops a well-defined limit-cycle state. The value of 
Fn drops below unity (i.e. becomes sub-Poissonian) for 
sufficiently strong coupling (k > 0.0075 in Fig. ^jp). 



C. Noise spectra of the system 

The steady-state of a system only gives information 
about average quantities. By also calculating noise spec- 
tra, further important information about the dynamics of 
the system can be obtained. We define the symmetrized 
noise spectrum for the two operators b and as, 



S,,,,iu) = lim / dr({6(t + r),6t(t)})e-- (7) 

J —OC 

= 23? lim / dT{{b{t + T),9{t)})e"^\{8) 



where b{t) = b{t) — (b), (b) = Tr[bpss] and 3? indicates 
the real part. The symmetrized noise spectrum has the 
property S'f, f,t(w) — Si,t i,{~uj). Note that we consider 
fluctuations about the steady-state of the system repre- 
sented by the limit f — > 00. 
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For system operators the correlation function can 
be evaluated by using the quantum regression theorem 
(QRT). The QRT states that the two-time correlation 
function can be written;^ 

lim (b{t + Tp{t)) = Trlbe^^b^pss] , r > 0. (9) 

In Liouville space this leads to,— 

S,,f>{uj) = m{{lo\Bn{uj)B^ro)), (10) 

where B is the Liouville space symmetrized superoperator 
whose relation to the Hilbert space operator b is given by, 



B\p{t)))^^{bp{t)+pm 



(11) 



and we define the pseudo-inverse of the Liouvillian, 

n{uj) = wi-iuj- cy^w, (12) 

with W an operator that projects away from the steady- 
state of the system, W = 1— |7'o))((/o|- The matrix formu- 
lation, given in Eq. (jlOp . is used to numerically evaluate 
the noise spectra in the next section. 

We can obtain helpful approximations and consider- 
able insight into how the resonator dynamics affect the 
noise by performing an eigenfunction expansion of the Li- 
ouvillian [Eq. ^] in TZ{u!) to obtain the alternative form 
for the spectrum, 



.p=i 



1 



-iuj — A, 



-{{lo\BK)){{l,\B^ro)) 



^ J. SR[Ap] jR[mP ] - (S[Ap] + c.)SK] 



5R[Ap]2 + (3[Ap]+u;)2 



(13) 



where = {{lo\B\rp)){{lp\B^ro)) . In almost all the cases 
considered here turns out to be real, and each term 
in the sum corresponds to a Lorentzian. However, in the 
case of the current noise we find that 3[»Tt^] 7^ 0, leading 
to a somewhat more complex lineshape)^ a point which 
we discuss further in Sec. El 

A knowledge of the eigenvalues of the system tells us 
where we can expect features in the spectra and how 
wide these features are. For the two separate systems of a 
damped oscillator and the SSET near the JQP resonance 
the eigenvalues are easily obtained. The eigenvalues for 
an oscillator coupled to a thermal bath are^^ 

X = -^mn-(l + lJm.\y,.., iTAlii, ±2, . . . (^4) 

The eigenvalues of the SSET do not have such a simple 
analytic form. For small Ej, the non-zero eigenvalues 
consist of a conjugate pair with a real part ~ — r/2 and 
i maginary pa rts close to the SSET frequency fissET = 
y^AE'^ + Ej/h, with the others real and of order — F. In 
the limit of zero coupling, the eigenvalues of the com- 
bined system will just be the sum and difference of these 
eigenvalues. 



III. SPECTRAL PROPERTIES OF THE 
RESONATOR 

The first resonator spectrum we consider is 

/oc 
dT{{a{t + T),aHt)})e'^- 
-00 

= m{{lo\ATZ{oj)A^\ro)), (15) 

where A\p{t))) 4=» ^ {ap{t) + p{t)a). The width of the 
peak that appears at; the resonator frequency in this spec- 
trum is the resonator linewidth, 70. For a superconduct- 
ing stripline resonator this spectrum can be inferred by 
probing the field that leaks out of the resonator— (e.g. 
via capacitive coupling to a transmission line^) , but here 
we will focus instead [in Sec.|V] on how the current noise 
spectrum can be used to obtain the resonator linewidth. 
For a resonator coupled only to a thermal bath this spec- 
trum consists of a single Lorentzian peak at the frequency 
of the resonator, with a width given by half the energy 
relaxation rate which in this case is simply 7ext/2. 

The behavior of S^ ,^i (lo) for AE detuned away from 
the center of the fc = — 1 resonance [Eq. ([5])] (where 
the resonator is in a fixed point state) and at the res- 
onance (with K such that the system is just inside the 
region where a well-defined limit-cycle exists) is shown 
in Figs. [3^ and [3Jd respectively. For the off-resonant 
case, peaks are observed at w = 0, a; ~ 17, w ~ f^sSET, 
w ~ — fissET and tu ~ —il + fissET- When at reso- 
nance we also observe the appearance of additional peaks 
at higher multiples of the resonator frequency. Also note 
that the peaks involving the SSET frequency have a much 
larger width since F ^ 7oxt ■ 

The spectra are dominated by the peak near to the 
frequency of the resonator, which is shown in detail in 
the insets. In each case the position of the peak is shifted 
slightly from the bare frequency of the resonator due to 
a renormalization arising from the interaction with the 
SSET. In both cases the peak remains Lorentzian. By 
comparing these two plots alone we see that the linewidth 
is much narrower for the resonant case. The spectra are 
shown for a relatively small value of k, but for increased 
coupling similar features are observed. 

We also calculate the spectrum of energy fluctuations 
in the resonator, 



Sn,nH = lim / dT {{n{t + T),n{t)}) e'^^ 



^m{{lo\Mn{uj)N\ro)). 



(16) 



where J^\p{t))) 4=> ^ {np{t) + p[t)n). For a resonator 
alone this spectrum consists of a single peak at w = 
with a width, 7„, given by the energy relaxation rate 
(7cxt in this case) and height 4(n^) /7cxt- For the cou- 
pled system the peak at w = is found to be Lorentzian 
to a very good approximation in both the limit-cycle and 
fixed point regimes, but not within the transition region. 
Wc show in Sec. IIVI that the width of this peak, 7„, is 



0.5 



1.5 2 2.5 



3.5 



FIG. 3: (Color online) Resonator spectrum 5„ j,t('^) for k = 
0.003 (other parameters are the same as in Fig. [2}. (a) off- 
resonance Ajf? = — 1.75eVds, f^sSET = l.lfi on a logj^Q scale, 
(b) on-resonance /S.E — — 1.59eVda, f^ssEX = f2 on a log^Q 
scale. The insets show the behavior around the peak at tj ~ O 
without logarithmic scaling (full lines) along with Lorentzian 
fits to the peaks (dashed lines). 



still given by the energy relaxation rate of the resonator 
in both the Umit-cycle and fixed point states. 

The behavior of the resonator Uncwidth. jq, and the 
rate 7„ obtained from Lorentzian fits to the appropriate 
spectral peaks are shown as a function of AE in Fig. U) 
The behavior far from resonance is well-understood: the 
resonator remains in the fixed point state and the SSET 
acts on it like an effective thermal batbi'^'^i^^ leading 
(for AE < 0) to a reduction in the energy relaxation 
rate and an increase in the effective frequency of the res- 
onator. Thus one expects that in this regime the res- 
onator linewidth should be equal to 7ri/2, as wc find. 
Close to the center of the resonance where the resonator 
is in a well-developed limit-cycle state there is a strong 
suppression of the linewidth, whilst the energy relaxation 
rate remains much larger. This behavior is precisely what 
one expects for a self-sustained oscillatori^ 

The fact that a given spectrum can be expressed in 
terms of the eigenfunction expansion of the Liouvillian 
as a sum of Lorentzians [Eq. [13] suggests that where the 
spectral peaks are found to be Lorentzian a single term 
in the expansion should dominate and the peak width 
be given by a single eigenvalue. Considering first the 
spectrum of energy fluctuations, S'„_„(w), we note that 
since this is the Fourier transform of the correlation func- 
tion {n{t + T)n{t)), we expect the places where (n^) and 
Sn,ni^) can be captured by a single term in the eigcn- 



FIG. 4: (Color online) Width of zero frequency peak in 
the Sn,n{i^) spectrum (upper curves) as determined from a 
Lorentzian fit (7n/2, solid line) and from the smallest non-zero 
eigenvalue (— Ai/2, dashed line). The lower curves show the 
linewidth of the resonator as determined from a Lorentzian fit 
to S'a at('^) (7f2, solid line) and from the eigenvalue (— 5RAn, 
dashed line), k — 0.003 with other parameters the same as in 
Fig. [2] Vertical dashed lines indicate the transition region. 



function expansion to be closely related. In the fixed 
point and limit-cycle states (but not the transition re- 
gion) the first term in the eigenfunction expansion de- 
scribes the variance in the resonator energy rather well^ 
(n^) ~ ((Zo|A/'|7'i))((Zi |A/'|ro)), where the corresponding 
eigenvalue Ai is the smallest non-zero one. The single 
term approximation to the spectrum is. 



5„,„(w)~43? 



1 



ioj — Ai 
4(-Ai)(n^) 



{{lo\U\nMimro)) 



(17) 



The validity of this expression is confirmed in Fig. 3] 
where we compare the eigenvalue Ai and 7„ and find 
excellent agreement except within the transition region, 
where the £7 = peak is not well-described by a 
Lorentzian and further terms in the eigenfunction expan- 
sion are required. 

For at it is expected that the eigenvalue closest to 
—ifl is the most important (this is certainly true for the 
decoupled system). We denote this eigenvalue Xn- For 
the eigenfunction expansion of the steady-state quantity 
({a, at}^ we find that (for fixed point and limit-cycle 
states and also in the transition region for weak coupling) 
the term corresponding to the Xq eigenvalue dominates, 
i.e. ^ + {n) ~ {{lo\A\rn)){{ln\A^\ro)). Thus wc expect that 
the spectrum around lu ~ Q should be well approximated 
by, 



Snai (w) 



4(-5ftAa) ((n) 



(18) 
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FIG. 5: (Color online) Widths of spectral peaks calculated 
from the eigenvalues as a function of k and /S.E. (a) logjQ [70], 
the logarithm of the linewidth of at('^) (b)7n, the width 
of the zero-frequency peak in Sn,n{ui), (other parameters and 
lines are as in Fig. [2} . 



where 51 is the renormalized frequency of the resonator, 
which is given by the imaginary part of the eigenvalue, 
Q.R = — 5[Ao]. The comparison of 70 with JRAn in Fig. 
13] shows convincing agreement. 

The calculation of the eigenvalues is much less numer- 
ically intensive than the calculation of the full spectrum, 
allowing us to calculate the linewidth for a large range 
of parameters. In Fig. [5^ we show the value of 70, as 
determined from the Aq eigenvalue, for the same range 
of parameters as Fig. [51 For weak coupling 7^ decreases 
with increasing (n) and hence the maximum in (n) cor- 
responds to a minimum in the linewidth. This behavior 
is analogous to that seen in a laser where the linewidth 
above threshold is given by the rate of phase diffusion;^ 
^laser _ ^ where G is the gain. However, a more com- 
plex behavior is apparent close to resonance for k > 0.005 
where (n) decreases with increasing k (see Fig. [2]) and 
the (single) maximum in (n) now corresponds to a local 
maximum in 70 with minima on either side. 

The behavior of manifests itself in the height of 
peak a,t uj ^ il in S'a^at(w), which from Eq. ^5)) . is 
~ 4 (n) /70 ((n) ^ 1 in the limit-cycle). The height 
of this peak as a function of AE initially shows a single 
maximum at AE = — 1.59eVds, which splits into a pair 
of maximas (with a minimum between) as the coupling is 
increased. Similar behavior was observed in the emission 
spectrum of a stripline resonator coupled to a SSET in 
a recent experiment.— However, this experiment was in a 
rather different parameter regime to the one considered 
here {Ej ^ HT) and in the experiment there were accom- 
panying features in the current, which we do not obtain 
here. A double peak in the spectral maximum as a func- 
tion of the detuning from resonance was also predicted 
in the emission spectrum for a closely related system in 



Ref. 0. In contrast to the linewidth, the energy relax- 
ation rate 7„ shown in Fig. [5]3 shows a similar behavior 
to that shown in Fig. [4] for all the values of k plotted. 



IV. ORIGIN OF THE PEAK WIDTHS 

We have shown that the widths of the peaks in the 
noise spectra are related to the eigenvalues of the Liouvil- 
lian. and that when the peaks are Lorentzian, the widths 
are given by a single eigenvalue corresponding to the 
dominant noise process. In this section we identify the 
eigenvalues with physical processes by comparing them 
to analytic expressions that describe those processes. 

We first consider 7„, the width of the uj = peak in 
Sn,ni^)- We do not expect the correlation function of n 
to depend on phase, and so in Section [TTI] we identified 
this term as the energy relaxation rate. This is supported 
by the fact that this term reduces to 7oxt in the zero- 
coupling limit. It is well known that when the resonator 
is in a fixed point state, the energy relaxation rate is given 
by the sum of the external damping and an additional 
effective damping due to the SSET^— and indeed we find 
that when the resonator is in a fixed point state, 7„ is 
well approximated by this total damping rate. We will 
now go on to show that 7„ also corresponds to the energy 
relaxation rate in the limit-cycle. 

In a limit-cycle, the effective damping can be general- 
ized to obtain an amplitude dependent damping^'^ '^^i'^^ 
7off(£') where E « (n) is the average energy. The cal- 
culation, which is described in Ref. [SJ, proceeds by de- 
riving a set of Langevin equations for the SSET and res- 
onator operators. Assuming then that the resonator en- 
ergy relaxes much more slowly than the SSET charge, 
the problem can be separated into two parts (though at 
the cost of neglecting some of the SSET-resonator cor- 
relations). Firstly the SSET charge dynamics can be 
solved treating the resonator as an harmonic drive with 
a fixed amplitude. The driven charge response as a func- 
tion of the resonator amplitude is then used to obtain 
an effective Langevin equation for the resonator alonei^ 
Writing down the corresponding equation for the aver- 
age amplitude of the resonator leads immediately to the 
effective damping, an expression for which (valid in the 
limit Ej <^ r) is given in the appendix. The average 
energy of the resonator obeys the equation of motion. 



— = -ME)E, (19) 

where jt{E) = 7cxt + 7eff(^') is the total damping rate. 
Although this equation is nonlinear, for small fluctua- 
tions about a stable limit-cycle with average energy Eq 
(given by 7t(£'o) = 0) we can linearize this equation to 
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obtain, 
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dE 
'dt 

7lin 



-/iin{E-Eo) 
dj{E)T 



dE 



(20) 



E — En 



where 7iin is the hnearized damping, which gives the en- 
ergy relaxation rate of the resonator near to a stable 
limit-cycle solution. The same approach also allows us 
to calculate the renormalized frequency of the resonator, 
rifl (as we discuss in the appendix). 

The calculation of jcsiE) require a^^i'^^ approximations 
in which certain correlations between the resonator and 
the SSET are dropped. Not surprisingly, this approach 
does not capture all of the coupled dynamics of the sys- 
tem and in particular it does not describe the small shift 
in the SSET frequency which arises due to the coupling 
to the resonator. This SSET frequency shift is apparent 
in the plot of (n) in Fig. [2^, where there is a clear change 
in the position of the peak as the coupling is increased 
(since this shift is much larger than than — 51, it can 
be attributed to a shift in ^sset)- The shifted frequency 
of the SSET can be obtained from a calculation which in- 
cludes more of the SSET-resonator correlations (for ex- 
ample from the shift in the SSET eigenvalues within the 
second order mean field equations described in Ref. 0) , or 
simply from the precise location of the peaks in the cur- 
rent noise spectrum of the SSET (Section |V|. In Fig. [6] 
we compare 7„ obtained numerically (from the Liouvil- 
lian eigenvalue Ai) to 7iin plotted against AE renormal- 
ized to take into account the shift in ft s set (for low- Ej 
and high-fi, the required renormalization is well approxi- 
mated by AE w AE + 2KeVds)- We find good agreement 
for these parameters which confirms that our identifica- 
tion of 7„ as the energy relaxation rate is indeed valid. 
The small difference between the curves in Fig. |6] arises 
because the expression for 7iin is strictly only valid in 
the small Ej limit; we find that a smaller value of Ej 
improves the agreement. 

We now consider 7n, the width of the u; = fl peak 
in S'(j_jjt(w). In the fixed point regime, we find that 
7n ~ 7ri/2, i.e. this width is also determined by the 
energy relaxation rate — indicating that the dominant 
fiuctuations of a are energy fluctuations in this regime. 
However, in the limit-cycle regime where the energy fiuc- 
tuations are small there is an additional process by which 
a can fluctuate, namely phase diffusion. The linewidth 
in a laser is determined by this phase diffusion rate^^i^i 
and we shall now show that this is also the case for the 
SSET-resonator system. 

Mapping the effective Langevin equation for the 
resonator— onto a c-number equation;^ we then follow 
the approach of Lax^i^ and derive effective equations 
of motion for the phase and amplitude. The calculation 
follows very closely that described in Ref. [H where the 
phase diffusion rate of a resonator coupled to a driven 
cavity is obtained. The effective amplitude and phase 




-1.6 -1.5 

FIG. 6: (Color online) Comparison of the width of the zero 
frequency peak in Sn,n{i^), as given by the eigenvalue 7„(Ai5) 
(solid line), with the energy relaxation rate 7iin(Ai5) (dashed 
line). Here k = 0.03, with other parameters the same as in 
Fig.H 



equations contain terms representing fluctuations aris- 
ing from both the environment and the the SSET island 
charge. These fluctuations (calculated in Ref. [H and 
summarized in the appendix) lead to effective diffusion 
terms for the amplitude {D^g) and phase (D^g) of the 
resonator about the limit-cycle, where the effective diffu- 
sion constants are valid on timescales long compared to 
the resonator period, but short compared to the damping 
rate. 

We find^ that as well as the direct phase diffusion 
we must also take into account another mechanism 
of phase diffusion. The frequency shift fl' = flu — de- 
pends on the amplitude, so any change in amplitude will 
lead to a change in frequencyi^ This means that a term 
describing amplitude fluctuations appears in the equa- 
tion of motion of phase, and thus there are both direct 
and indirect contributions to the phase diffusion. 

We combine the terms arising from direct phase dif- 
fusion and the contribution from amplitude fiuctua- 
tions (neglecting cross-correlations) to obtain the total 
linewidth 7^ = 7^ + 70> 



14' 



D,^, + D+,{E) 



4(n) 



7lin 



D,^t + D-g{E) 



(21) 



dn'(E) 
dE 



is the frequency shift lin- 



at zero 



where = i?o 

earized about the limit-cycle, and -Doxt 
temperature. In Fig. [7] we compare 7^ with the numer- 
ically obtained result for the linewidth from the eigen- 
value, 7n. We see that near to the center of the reso- 
nance the indirect diffusion term tends to zero, but it 
contributes significantly when the system is not exactly 
on resonance. 

Looking at Fig. [7] we see that the double-peak structure 
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-1.8 -1.7 -1.6 -1.5 -1.4 -1.3 



FIG. 7: (Color online) Comparison of the linewidth 7n of 
at {^) obtained from the An eigenvalue (plotted as a func- 
tion of Ai5) with the analytic expression 7^ = 7^+7^ (plotted 
as a function of A.E). Also shown are the direct (7^) and in- 
direct (7J) phase diffusion contributions. Here k = 0.03, with 
other parameters the same as in Fig. [S] 

seen in 7^ at strong couplings arises from the combination 
of the maximum in 7^ at the center of the resonance 
and the strong increase in 7^ off-resonance. At weaker 
couplings the limit-cycle regime becomes much narrower 
as a function of AE. Here we find that the combination 
of a less pronounced maximum in 7^ and a much sharper 
minimum in 7JJ together lead to a single minimum in 7^. 

V. CURRENT NOISE 

In this section wc investigate the extent to which fea- 
tures in the resonator noise spectra from Sec. IIIII mani- 
fest themselves in the current noise. The low frequency 
limit of the current noise in this system was discussed in 
Ref.ll, but here our interest is in the dynamical timescales 
which become imprinted in the widths of the peaks in 
the finite frequency current noise. The finite frequency 
current noise through the SSET can be split into contri- 
butions from the current noise at the two junctions and 
the charge noise of the island,— 

SiA^) = lsi,j,{u;) + ^Si,j^ico)^\Lu^SQ^Q{Lo), (22) 

where we have assumed equal junction capacitances. We 
choose the Cooper pair tunneling to take place at the 
left hand junction. To calculate the charge noise and the 
current noise at the left hand junction we can again use 
the quantum regression theorem as described in Sec. Ill C[ 

Si^j^iio) = m{{lo\lLn^)lL\ra)), 
SQ,Qiio) = 43f^((/o|Q7^(a;)Q|ro)), (23) 



where the current and charge operators arc defined as, 

I^=^e^{\2){0\^\0){2\), 
Q = e(|l)(l|+2|2)(2|). (24) 

and the superoperator forms are defined in the same man- 
ner as Eq. ([TT]). For the current noise at the right hand 
junction a counting variable approach can be usedi^ The 
resulting expression is very much the same but with the 
addition of a self-correlation term, 

Si^j,{lj) = 2e{{lo\lR\ro))+m{{lo\lRn{Lu)lR\ro)), (25) 

where the current super operator for the right hand junc- 
tion is defined, 

Inlpim = er(|0)(l|p(t)|l)(0| + |l)(2|p(i)|2)(l|). (26) 

Wc give results in terms of the current Fano factor, 
Fi{uj) = Sij{uj)/2e{I). For the SSET alone the cur- 
rent noise spectrum, for these parameters, consists of a 
large peak at w ~ f^ssET and a much smaller peak at 
w = 0. 

In Figs. [5^ and[5jD we show the current noise spectra 
for the same off-resonant and resonant parameters used 
in Figs. [3^ and[3)D. The SSET current has a non-linear 
dependence on the resonator position^ leading to peaks at 
higher harmonics?*^ oi fl. In the off-resonant regime the 
resonator state consists of Gaussian fluctuations about 
a fixed point and only the tu ~ 217 harmonic is seen in 
the current noise^i^ and it is very much weaker than 
the peak at w ~ fi. On-resonance, the resonator motion 
consists of oscillations of a (relatively) much larger am- 
plitude and hence leads to clearly visible peaks at higher 
harmonics of fl. 

It is interesting to note that for the peak near cj ~ 
the off-resonant current noise spectrum displays an im- 
portant difference from the resonator spectra: as is shown 
in the inset (Fig. [S^) the peak is not of a Lorentzian form. 
The non-Lorentzian shape of this peak in the fixed point 
regime can be understood^ as arising from an extra el- 
ement of correlation between the SSET and resonator: 
the SSET current picks up the fluctuations in the res- 
onator motion which were originally driven by the SSET 
charge motion (rather than a totally uncorrclated exter- 
nal bath)i^ This effect is rapidly reduced with increas- 
ing external damping and temperature. Note that in the 
limit-cycle (resonant) region the peak at ~ 12 is again 
a Lorentzian. 

Nevertheless, we find that the non-Lorentzian peak at 
w ~ J7 in the fixed point regime can still be fitted by a 
term in the eigenfunction expansion Eq. [T3] (but with a 
complex matrix element m^), and hence the parameter 
describing the width, 3?An,/, can be obtained over the 
whole range of AE. Comparing 'StXaj obtained in this 
way with the linewidth 7^, and comparing the width of 
the uj = peak in Sijiuj) with 7„ for the parameters 
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FIG. 8: (Color online) Finite frequency current noise spec- 
trum for K = 0.003 with the other parameters the same as in 
Fig-E] (a) Off-resonance Ai5 = — 1.75eVds, (b) on-resonance 
for /S.E = — 1.59eVds- The insets show expansions of the nar- 
row peaks at ~ f2 without logarithmic scaling (full lines) 
along with fits to the peaks (dashed lines) with the form of 
Eg. 1131 inset (a), and Lorentzian form, inset (b). 

in Fig. |4] wc find in both cases that the pairs of curves 
overlay one another (the resulting plot thus effectively 
duplicates that shown in Fig. 3] and hence is not shown 
here) . Hence the current noise provides a direct measure 
of the spectral properties of the resonator spectra and, 
in particular, gives access to the energy relaxation and 
phase diffusion rates. 



VI. CONCLUSIONS 

We have studied the noise properties of a resonator 
coupled to a SSET in the vicinity of the JQP resonance, 
focussing our analysis on the case of a resonator whose 
frequency is large compared to the quasiparticle decay 
rate. We analyze the resonator spectrum that includes 
fluctuations in both the resonator energy and phase, 
Saa^{^); "^c\\ as the spectrum of energy fluctuations 
S'„ „t(w). The main feature in the spectrum S'(j.at(w) is 
a large peak at the resonator frequency, w while 
„t (w) is dominated by a peak at w = 0. We find that 
outside transition regions, these peaks are Lorentzian 
with widths, 70 (the linewidth) and 7„, respectively, 
which are each controlled by a single eigenvalue of the 
Liouvillian. Examining the behavior of these eigenvalues 
for a range of coupling strengths, we find that for large 
enough coupling Aq has a double minimum as a func- 



tion of the detuning from resonance. An analysis of the 
coupled dynamics of the system allows us to identify the 
peak width 7„ with the energy relaxation rate and 70 
(within the limit-cycle regime) with the phase diffusion 
rate. We show that the double minimum in the linewidth 
arises from a combination of direct and indirect phase dif- 
fusion. Finally, we show that these rates can be extracted 
from the current noise, and thus the current noise pro- 
vides a direct measure of the key dynamical time scales 
of the resonator. 



Acknowledgements 

This work was supported by EPSRC under Grants 
EP/D066417/1 and EP/C540182/1. 



APPENDIX A: ENERGY DEPENDENT 
DAMPING AND DIFFUSION 

In this appendix we present further details on the ana- 
lytical approximations used to obtain the effective damp- 
ing and diffusion rates. The details of the methods used 
in the calculations (and the approximations required) 
can be found in Refs. [13 and fSQ. Although the expres- 
sions obtained are only valid for small Josephson energy, 
Ej ^ hT^ they are valid for all values of O/F, i.e. for both 
fast and slow resonators. The basis of the calculation is 
that the resonator energy is assumed to relax very slowly 
compared to both the resonator period and the quasi- 
particle tunneling rate, so that its effect on the SSET 
charge can be approximated as a harmonic drive. The 
behavior of the SSET under the influence of this drive 
is then found, and the result then fed back to obtain a 
coarse-grained equation of motion for the resonator. 

The effective damping, 7off, and frequency shift, fi' = 
VLji — 57, for a resonator with average (real) amplitude 
\/E are given by,— 

l.s{E) . x,n^^E]V^ (3r + 2z») 

(Al) 

where Xq = y^h/{2mfl), the Bessel function of the flrst 
kind Jn{z) is a function of the scaled amplitude z = 
2\fExsl Xq and we define = ^ + ?'(Jln + /S.E/h). 

In a similar way, the fluctuations of the charge on the 
SSET— lead to effective diffusion constants for the am- 
plitude and phase of the resonatori^ The details of the 
calculation are rather involved, and the calculation fol- 
lows the same approach as that given in Ref. [s^ hence 
here we merely summarize the results. The total diffusion 

rate is given by D^^ = -^^tt^VIk'& {of + + Of), 
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where £>J^ refers to the phase diffusion, and the D^f^ 
terms to the ampHtude diffusion. The first term is, 

1 V- |J„+iWCi± J„_i(z)Ci*|' 



where 



3r + 2in 



The second term is given by, 

j^±_Y- {Jn+l{z)Ci ± Jn^l{z)Ct) 
2 / ^ h* 



^ ' J„+i(z)C| ^ J„-l(z)C2 



h* 



h* 
"■n-1 



with 



Co 



2r + ii7 
(r + if])2 



The final term is given by. 



hn 



± C4J„-l(-2)Jn+l(2) 



'n+1 



where 



C4 



5r2 + 2f]2 



T + i2n (r + jsi)2(r + i2r2)2 



The effective damping and diffusion can now be used to 
calculate the linewidth. These quantities can in addition 
be integrated to give an effective potential, allowing a 
calculation of the resonator distribution.— We note that 
this result generalizes the phase diffusion rate calculated 
m to the regime where il/T is no longer small. 
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